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, We present a density functional theory (DFT) calculation in the generalized gradient 

■ approximation to study the possibility for the existence of Jahn- Teller (JT) or pseudo 

' Jahn- Teller (PJT) type local distortions in the La-Ba-Cu-O superconducting system, 

^—j- ' We performed the calculation and correspondingly group theory classification of the 

£f*) , electronic ground state of the Cu06 — 10 elongated octahedra cluster, immersed in a 

f - ^ , background simulating the superconductor. Part of the motivation to do this study is 

' that the origin of the apical deformation of the CuOg - 10 cluster is not due to a pure JT 

f ^ ' effect, having therefore a non a priori condition to remove the degeneracy of the electronic 

ground state of the parent regular octahedron. We present a comparative analysis of the 
symmetry classified electron spectrum with previously reported results using unrestricted 
Hartree-Fock calculations (UHF). Both the DFT and UHF calculations produced a non 
degenerate electronic ground state, not having therefore the necessary condition for a 
pure JT effect. However, the appearance of a degenerate E 9 state near to the highest 
occupied molecular orbital in the DFT calculation, suggests the possibility for a PJT 



effect responsible for a local distortion of the oxidized CuO g 9 cluster. 
Keywords: Cuprates superconductors; density functional theory; Jahn- Teller effect 



1. Introduction 

Since the discovery of high T c superconductivity, the understanding of the mecha- 
nism responsible for this novel phenomenon has been one of the most challenging 
problems in condensed matter physics. To date, it is known to involve the pairing 
of charge carriers, although the precise nature of these carriers and the mechanism 
by which the pairing occurs is still an open problem. 1 ~ 4 The evidence that strong 
clcctron-phonon interaction in oxides as well as in mixed-valence systems can occur 
owing to polaron formation, together with the fact that a possible mechanism for 
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polaron formation is the Jahn- Teller (JT) effect, 5 where spontaneous lattice dis- 
tortions remove the degeneracy of the electronic ground state, led to Bednorz and 
Miiller 6 to the discovery of high temperature superconductivity. One reason to sus- 
pect the JT effect might play an important role in those systems, is the fact that 
all superconducting copper oxides have perovskite-like crystal lattices, with Cu 2+ 
ions surrounded by 2 ~ ions. One of the most studied system, La-Ba-Cu-0 has 
this particular structure, containing sheets of corner-sharing CuOg octahcdra. The 
Cu-0 planes, where it is believed that superconductivity and charge transport are 
mostly confined, have an electronic structure, within the ionic model, built up from 
the 3d 9 states of Cu 2+ and the 2p 6 states of 2 ~. 

In the insulating parent compound of the La based superconductor La2Cu04, the 
CuOe octahedra are elongated, with a Cu-0 distance of 1.90 A within the plane and 
2.42 A perpendicular to it. This distortion, some authors believed, could be at least 
partly due to a static JT effect. 7 This can be seen from the fact that for the regular 
octahedron we have the Oh symmetry, with a splitting of the d 9 electron states into 
E g + T 2fl irreducible representation (IR). The corresponding decomposition of the 
ions vibrational modes Qk, according to the IR of the O^ group gives T — &i g + 
c g+ t2 S + 2ti u + t 2 „. The general condition in order to have a JT effect, is that the 
generalized force between nuclei, defined as the matrix element of the force operator 
(-dV/dQk), between the total electronic ground state fa, transforming respectively 
as the Tk and Ti IR of the corresponding symmetry group, should be different from 
zero. From general group theory rules, we know that this condition can be fulfilled 
if r* x T k x Tj C Ti is satisfied (or equivalently T* xT, C r fe ), where Ti is the 
identity IR. Hence for the regular octahedron having the Oh symmetry group, the 
JT condition is fulfill, being the e g and t2 9 modes the possible JT active modes. 
Hence in order to have a pure JT effect and as a consequence a completely removal 
of the electronic degeneracy, the Cu site symmetry should be cubic. This would be 
precisely the case for an isolated regular CuOg octahedra. However in the insulating 
parent compound of the La based superconductor La2Cu04, the Cu site symmetry 
is far from being cubic. Even if the CuOe octahedra would be regular, the crystal 
field would not possess cubic symmetry. It has also been shown 8 (and references 
cited there), that at least half of the octahedral distortion in the CuOe cluster 
can be accounted for, as a consequence of the strongly layered crystal structure 
and the large ionic interaction. From the above mentioned arguments one would 
expect that the partially filled electronic E ff state in the regular CuOe octahedra 
might only be partly responsible for the CuOe elongated shape, having therefore 
a non a priori condition for the complete removal of the electronic degeneracy. As 
a clear example of an apical distortion where the JT effect does not play any role 
is when the Cu ion is substituted by Ni, obtaining La2Ni04. This system has the 
same structure as La2Cu04, but the Ni+ 2 is not a JT ion. There are also in the 
literature some proposed models where the JT effect plays a definite role in high 
T c superconductivity beyond that of producing the static deformation of the CuOe 
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octahedra. 9-11 

Experimental results in Cu K-edge absorption spectroscopy XANES (x-ray 
absorption near edge structure) and EXAFS 12 (extended x-ray absorption fine 
structure) as well as PDF (pair-distribution-function) analysis of neutron powder- 
diffraction data 13 have shown the existence of important inhomogeneities in 
the high-Tc superconductor structures. In particular, the CuOe octahedra in 
Lai.85Sro.i5Cu04 exhibits local deformations in the apical direction as well as in the 
planar Cu-0 bond length due presumably to changes in the local charge state. 12 
This fact together with the previously mentioned arguments concerning the pos- 
sible symmetry of the electronic ground state, raises the question as whether this 
deformation is somewhere connected to a JT or pseudo Jahn- Teller (PJT) kind 
of distortion. In order to determine whether the JT or PJT effect might play any 
significant role on the local distortion on the La-Ba-Cu-0 system, it is necessary 
to determine the symmetry properties of the total electronic ground state. Being 
the condition for a JT that this state must be degenerate. If it is not the case then 
the possibility of a PJT effect emerges and it is determined by no null matrix ele- 
ments of the force operator (-dV/dQk), between the non degenerate ground state 
and closed excited states. Concerning the PJT effect it has been pointed out 14 that 
there are no a priori exceptions for vibronic coupling effects. The presence of elec- 
tronic degeneracy is just an important special case of vibronic coupling. Bersuker 15 
has introduced a PJT vibronic coupling through the mixing with the 1 Bi g excited 
state to explain a polaron formation mechanism in high T c superconductors. It is 
worth mention that the PJT effect has been some times misinterpreted as a second 
order correction, and in fact this effect may be very strong, specially when the JT 
effect is absent. The study of those effects might be relevant in understanding the 
electron-lattice interaction. 

In this paper, we examine the possibility of having a local distortion in the 
CuC>6 octahedra in La-Ba-Cu-0 superconducting system originated by a JT or 
PJT effect. This is done under the assumption of a Cu06~ 10 cluster model for 
the superconductor. The CuC>6~ 10 cluster is immersed in a set of static charges, 
located in the proper positions in order to simulate the crystalline background. 
This background stabilizes the cluster and preserves the point symmetry group. 
The electronic spectrum is calculated using the density functional theory in the 
generalized gradient approximation (DFT-GGA), and classified according to the IR 
of the symmetry group D^. 

The paper is organized as follows. In Sec. 2 we introduce the cluster model 
representing the La-Ba-Cu-0 system, and the methods used for the calculation 
and symmetry classification of the electronic spectrum. In Sec. 3 we present the 
results and discussion for the spectrum obtained and compared with the results of 
a unrestricted Hartree-Fock (UHF) calculation, and discussed the possible JT or 
PJT type of distortions. The conclusions are presented in Sec. 4. 
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2. Cluster Model 

Embedded cluster models are used to describe real systems by treating the cluster 
electronic structure with some degree of sophistication by using ab initio or semi- 
empirical methods, where as cruder approximations have to be made in order to 
treat the background used to embed the cluster. The basic philosophy in the use of 
cluster models is to study essential characteristics of the solid which are primarily 
determined by local properties that can be simulated by cluster calculations. Its use 
has been mainly justified by many successful applications including, by example, 
the La2CuC>4 undoped cuprate 16 ' 17 or the L^-xSr^CuC^ superconductor system. 18 
In this work we used the cluster model to study the effect of the localized orbital 
electronic states on the possibility for the existence of Jahn- Teller (JT) or pseudo 
Jahn- Teller (PJT) type local distortions in the La-Ba-Cu-0 superconducting sys- 
tem. 

We study the elongated octahedral CuC>6 -10 cluster immersed in a point 
charges environment according to the structure of the tetragonal superconduc- 
tor (La,Ba)2CuC>4 with the experimental lattice constants: a = 3.7873 A and 
c = 13.2883 A and with a Cu-0 apex distance of 2.42 A. 19 This background 
simulates the crystalline surrounding. 

The necessity to preserve the experimental bond lengths of the cluster, lead 
us to adjust the background arrangement in order to stabilize the cluster system. 
The values for the point charges environment were determined by performing a 
UHF cluster stability analysis. 20 The first trial of stability calculations was done 
surrounding the cluster with fixed point charges corresponding to the formal charges 
values La 3+ , Cu 2+ , 2 ~, and placing 1 1538 ~ at the outer boundary (in a spherical 
radius of 8 A ) accumulating 178 charges which neutralize the -10 charge of the 
cluster. We found this structure unstable in this background. However, if we change 
the formal charge of the nearest copper atoms in Cu-0 planes to Cu 1+ (value 
consistent with possible Cu valence states), we recuperate stability with an absolute 
value for the energy gradient of 2.2 xl0~ 3 a.u. for the oxygen atoms in the plane. 
At this stage we did not try to make a full structure optimized calculation. 

Recent works dealing with similar cuprate clusters have shown that the use of 
pseudopotential models for the nearest atoms to the cluster are advantageous in or- 
der to prevent an artificial polarization of the oxygen electrons towards the (in other 
way) bare positive charges. 17 ' 18,21 In a previously reported work, 20 we performed 
an electronic structure calculation for the neutral CuOg^ 10 cluster embedded in 
the modified point charges environment using different basis in the HF formalism. 
Looking at the hybridization of the different wave functions, we do not observe the 
artificial polarization effect. This is due to the fact that we are using an effective 
charge of +1 for the nearest Cu atoms on the Cu-0 planes . The only case where 
we note anomalies was in the use of minimal gaussian basis 3-21G, where we got a 
wrong symmetry for the total wave function. By the addition of a polarization term 
to the oxygen 3-21G basis, we recuperated the correct wave function symmetry. It 
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should be mentioned that the point charges substitution by the pseudopotentials 
does not guarantee the stabilization of the structure without a full relaxation of the 
cluster. In a DFT structure optimization calculation using bare pseudopotentials, 
relaxing only the apical length, Plibersek et al. 22 found a +0.04 A shift from the 
correct experimental value for the apical distance Cu-O. Hence, since our main in- 
terest is in the analysis of the symmetry properties of the electron spectrum we try 
as a first model calculation a background with point charges. 

When a point charges set is used to simulate the crystalline environment, some 
small adjustments in the charge values or its positions, or the inclusion of a ghost 
charges set are made in order to obtain the Madelung potential on the points corre- 
sponding to the atoms positions in the cluster. We do not assume this procedure here 
because the difference between the potential resultant from the number of charges 
in our model and the Madelung potential becomes rather small. 23 Even more, as- 
suming this procedure, the background symmetry could be changed if some charges 
are moved from its original positions. In our model we need to preserve the point 
symmetry of the environment in order to obtain the correct IR's corresponding to 
the cluster molecular orbital (MO) and the nature of the vibronic coupling. 

In summary, our model consists of a Cu06~ 10 cluster embedded in an environ- 
ment formed by point charges placed at the positions of the nearest shell of atoms 
outside the cluster with the outer shell always being oxygen and the immediate cop- 
per atoms in the Cu-0 plane with a charge +1. For the present calculation we have 
surrounded the cluster with 178 fixed point charges corresponding to La 3+ , Cu 2+ , 
2 ~, Cu 1+ , and O 1 026 ~ at the outer boundary, with a net charge of +10 which 
neutralized the -10 charge of the cluster. The radius of the outer shell was fixed at 
8 A . In order to consider the possible undesirable basis effects on the electronic 
spectrum we selected a 4-2 1G basis with polarization d-components in the oxygens 
and no diffuse components s or p were added. The additional oxygen d-components 
privileges the Cu3d-02p hybridization over the artificial polarization 02p-(Cu 1+ , 
La 3+ bare charges). The arrangement selected, guarantees approximately the sta- 
bility of the cluster, as well as the symmetry properties of the electronic states 
and small polarization effects. 

In a simple ionic model, the molecular levels formed by the 3d 9 states of Cu 2+ 
hybridized with the 2p 6 states of 2 ~ immersed in a D 4 /j crystalline field pro- 
duce a nondegenerate electronic ground state. The same result is obtained within a 
quasi-molecular approximation in an extended Huckcl model and also in open shell 
restricted and unrestricted Hartree-Fock (RHF-OS and UHF) calculations for the 
electronic ground state. 20 This shows that within these approximations we do not 
find a possibility for pure JT distortions. The question is whether the same sym- 
metry for the total electronic wave function is preserved in an electronic spectra 
obtained from ab initio calculations where higher electron correlations are taken 
into account. In order to investigate this situation we have performed a DFT-GGA 
calculation to determine the electronic ground state for the apical octahedron cluster 
together with their symmetry classification. This allows us to analyze the possibility 
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of realization of total degenerate electronic wave function. 

The DFT calculation for the electronic ground state was performed us- 
ing The Cambridge Analytic Derivatives Package CADPAC 6.1 24 running in a 
CRAY/YMP-4. This calculation was done using the generalized gradient approx- 
imation (GGA) 25 ~ 28 with a 3-21G basis set for Cu and a 4-31G* for O (where 
* means polarized). The GGA improves the previous commonly used DFT-LDA 
(local density approximation) incorporating some nonlocal effects to the exchange- 
correlation functional using the gradient of the electron density. We employed a 
convergence critcrium of lxlO -7 for the density. The expected value for the spin 
operator S 2 in the ground state, was systematically checked to be 3/4. The sym- 
metry of the total electronic wave function and hence of the density is guaranteed 
by the symmetry of the hamiltonian (D^). For the exchange-correlation X-C func- 
tionals, the BP86 approximation (Becke 29 exchange and Perdew correlation 27 ) was 
used. We consider a system formed by 87 electrons of the CuOg -10 octahedron 
cluster which corresponds to the Cu 2+ and 60 2 ~ ions. For the calculation, a total 
of 119 Gaussian base functions were employed, corresponding to 29 for the Cu 2+ 
and 15 for each of the six 2 ~.The number of a and (3 occupied states are 44 and 
43 respectively, corresponding to the 87 electrons within the cluster. The obtained 
electronic wave functions were classified according to the IR of the symmetry 
group using a simple algorithm based on the projection of characters that permits 
to find the classification according to the IR of a given symmetry group with the 
application of at most three symmetry operation. 30 

3. Results and Discussion 

We have studied within density functional theory in the generalized gradient approx- 
imation, the electronic ground state for the CuOg -10 elongated octahedra cluster. 
This cluster represent one of the T c superconductors that have been suggested as a 
prominent candidate to show the JT effect : La-Ba-Cu-O. 

The electronic spectrum for a states in the DFT for the Cu-3-21G plus O- 
4-21G* basis set in joint to the UHF spectra for Cu-3-21G plus 0-3-21G*, 0-4- 
21G* and 0-6-31G* is shown in Fig. 1. The HOMO (highest occupied molecular 
orbital) states correspond, in both UHF and DFT, to a nondegenerate a level 
whose IR is 4i?i 9 giving, after making the product with the (3 electronic wave 
function, a total orbital nondegenerate electronic ground state wave function whose 
symmetry is 2 B\ g . It should be mentioned that in Ref. 21, although the authors do 
not explicitly refer to the symmetry of the ground state, it is possible to infer it from 
their components of the atomic orbitals, corresponding the 2 Bi ff , which agrees with 
our result. As can be seen from the spectra, all the HOMO states correspond to 
nondegenerate B\ g states. Also all a-LUMO (lowest unoccupied molecular orbital) 
states are the nondegenerate IR HA\ g . It should be noted that the a-(HOMO- 
LUMO) gap is reduced considerably in DFT respect UHF, from 11.02 cV to 2.4 eV 
for the same basis. This result is consistent with the well known fact that the HF 
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Fig. 1. DFT electronic spectrum for the a states using the Cu-3-21G plus 0-4-21G* basis set. 
This is compared with the UHF spectra using the Cu-3-21G plus 0-3-21G*, 0-4-21G* and 0-6- 
31G* basis sets respectively of the CuOg - 10 cluster. 

method overestimates the HOMO-LUMO gap in many molecular systems. 

In Fig. 2 we show the DFT spectra for a and (3 states and in Table 1 the 
contributions from the copper 3d and the planar and apical oxygens 2p to the 
highest occupied KS-BP86 MO's. The LUMO is the f3 state 4B lg located 1.8 eV 
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Fig. 2. Spin-unrestricted DFT-GGA spectrum for the Cu06 10 cluster using the Cu-3-21G plus 
0-4-21G* basis set. 

over the HOMO, compared with 1.4 eV reported in Ref. 21. We would like to 
remark some of the particular features shown in the DFT spectrum. There are two 
molecular orbital "bands" in the occupied electronic spectrum under the HOMO, 



February 2, 2008 1:8 WSPC/INSTRUCTION FILE JTdflJMP 



DF electronic spectrum of the CuOq 10 and JT distortions in La-Ba-Cu-O. 9 



Table 1. Square of the copper 3d and the oxygen 2p (planar and apical) expansion coefficients of the highest 
occupied KS-BP86 MO's of the Cu0 6 ~ 10 cluster. 
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separated by a 2.9 eV gap. The first band with energy between and -0.91 eV 
have MO's with predominantly 3d cooper components: 4B lg (HOMO), 3E ff , 2B 2g , 
and 10Ai g . The second band with energy between -3.81 and -5.10 eV include MO's 
composed mainly by oxygen 2p components: lA 2ff , 7E M , 1B 2m , and 2E g (Table 1). 
This feature is not present in the UHF spectrum, where this MO bands are mixed. 
Similar results were obtained in Ref. 21: the first band between and -1.5 eV and 
the second band between -2.6 and -4.4 eV. 

We show in Fig. 3 the predominant d orbitals for the 4B iff (HOMO), 3E g , 2B 2g , 
and 10Ai s containing the 3d 9 states of Cu 2+ . In the HOMO state 4Bi s , the 3d :E 2_j / 2 
state of Cu is hybridized with the 2ptr O orbitals along the Cu-0 bond in the x-y 
plane. This state is occupied only by one a electron and is orbitally nondegenerate. 
It should be noted however that it is spin degenerate, but since it has a degeneracy 
of order two in spin, it does not produce a JT deformation. 31 It is basically the 
same kind of hybridization we got for the UHF wave function, except by a lower 
contribution from the 2ptr O orbitals. The occupied state nearest to the HOMO, 
the double degenerate orbital 3E ff , could be seen as the result of hybridization of 
the 3d xz {3d yz ) states with the p x and — p^ (p y and — p y ) oxygen orbitals along the 
apical axis with no contribution of the p z and -p z oxygen orbitals along the x (y) 
axis. The 2B 2ff orbital becomes the hybridization of the Cu 3d xy state with the p7r 
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Fig. 3. Hybridzation of the DFT-GGA wave function for the predominantly d orbitals of the 
Cu0 6 ~ 10 cluster: (a) AB lg (HOMO), (b) 3E g , (c) 2B 2g , and (d) WA lg . 

orbitals of the oxygens in the x-y plane (perpendicular to the Cu-0 bond). Finally, 
the 10Ai g orbital is the hybridization of the 3d z 2_ r 2 state of Cu with the 2pcr states 
of the oxygens in the apical z axis and with the p x , —p x , Py and — p y orbitals of 
the oxygens in the x-y plane along the Cu-0 bonds. 

Hence, since the electronic ground state is the nondegenerate 2 Bi 9 state, we can 
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not have a residual static JT effect in the Cu0 6 ~ 10 cluster. The possibility for a PJT 
deformation, in this state of oxidation, is low since the HOMO-LUMO gap which 
is about 2 eV is much greater than the characteristic vibrational energy (20-100 
meV). However, in the oxidized CuOg -9 cluster, the spin-singlet 1 Ai g electronic 
ground state could be mixed with some near excited states giving non-zero PJT 
vibronic couplings. Those states could be the spin-singlet (predominantly copper 
3d states) and/or 1 Bi ff . This result is consistent with the assumption made 
in the model introduced by Bersuker 15 to explain a polaron formation mechanism 
in high T c superconductors. There he assumes a PJT vibronic coupling through 
the mixing with the 1 Bi g excited state. Here we have shown that besides the 1 Bi g 
excited state, the inclusion of the 1 E g excited state in a PJT vibronic coupling might 
be also relevant. In the context of the model we have presented in this work, the 
PJT effect may produce a coupling with the e g and/or bi s vibrational modes which 
could be responsible for the local deformations in the superconducting system. 

Recent DFT Becke-3-Lee- Yang-Parr band calculations 32 in the I^-nSr^CuC^ 
doped system support a highly inhomogeneous hole formation with symmetry 1 Ai ff 
adjacent to the Sr impurity, giving us additional reliability about our final conclu- 
sion. 

4. Conclusions 

We have performed ab initio calculations in the CuOg -10 elongated octahedra clus- 
ter having the symmetry group D^, representing the La-Ba-Cu-0 system, to in- 
vestigate the possibility of having a Jahn- Teller or pseudo Jahn- Teller effect as the 
origin for the local distortions found in this system. This study has been motivated 
by the fact that the apical distortion of the CuOg octahedra in the insulating par- 
ent compound of the La based superconductor La2CuC>4 is only partially due to a 
static Jahn- Teller effect, having therefore a non a priori condition to remove all the 
degeneracy of the electronic ground state of the parent regular octahedron with Oh 
symmetry. 

The DFT-GGA electronic spectrum was compared with the corresponding UHF 
calculated in a previous work. As a result in both cases, we found that the electronic 
degenerate states correspond only to intermediate levels, which are fully occupied. 
The symmetry of the total electronic ground state wave functions being therefore a 
nondegenerate 2 Bi g state, ruling out the possibility of a residual static JT effect in 
the CuG-6 -10 cluster. It should be noted however, that the DFT electronic spectrum 
contains a degenerate 3E S state very close to the HOMO 4Bi 9 state. This belongs 
to a first set of filled molecular levels, including the HOMO state, which are spread 
on a width of the order of 1 eV, and which are characterized by being predominately 
formed by the 3d states of Cu. This set is separated by a gap of the order of 3 eV 
from another set of molecular states originated mainly by 2p states of O. 

Since the HOMO-LUMO gap (about 2 eV) is much greater than the character- 
istic vibrational energy (20-100 meV), the possibility for a PJT deformation, in this 
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state of oxidation is low. However, in the oxidized Cu0 6 ~ 9 cluster, the spin-singlet 
1 Ai g electronic ground state could be mixed with some near excited states giving 
non-zero PJT vibronic couplings. Those states could be the spin-singlet (predom- 
inantly copper 3d states) 1 E ff and/or 1 Bi 5 . We therefore conclude that the PJT 
effect in the Cu06 -9 cluster may produce a coupling with the e 9 and/or bi g vi- 
brational modes. This PJT effect could be the responsible for the local distortions 
observed in the CuC>6 octahcdra in the Lai.ssSro.isCuC^ superconductor. 
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